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FIG. 1 Creation of an entangled electron-hole pair by appli- 
cation of a voltage difference V between two metals separated 
by a tunnel barrier. The Fermi sea consists of the filled states 
in the conduction band. The spins of the electron (e) that 
has tunneled and the hole (h) it leaves behind are entangled 
in the state 2" 1/2 (| ThTe) + I UU))- 



in the conduction band of a metal are filled with elec- 
trons up to a maximal energy. This "sea" of electrons is 
called the Fermi sea and the maximal energy is the Fermi 
energy. 1 An unfilled state in the Fermi sea behaves like 
a positively-charged electron, called a hole. 

Electron-hole excitations in the Fermi sea are created 
by application of a voltage V over an insulating barrier, 
as illustrated in Fig. ^ An electron can tunnel through 
the barrier, leaving behind a hole. The Fermi energy Ep 
is higher by an amount eV on one side of the barrier 
-than on the other, so that the hole is elevated to the 
■fame energy as the electron. 

A few years ago it occurred to us that the spatially 
separated electron-hole pair created by a tunneling event 
is a spin- entangled Bell pair This simple observa- 
tion, that single-particle tunneling in a metal produces 
entanglement, came as a surprise. Earlier proposals for 
cntanglcrs of spatially separated electrons (reviewed in 
Refs. 0,11) had relied 

on electron-electron interactions 
in one way or the other: the Coulomb interaction in a 
quantum dot the (phonon- mediated) pairing in- 

teraction in a superconductor 0, 0, [nj, 0, ^ , or the 



A. Preface 



Enrico Fermi, to whom this School of Physics is dedi- 
cated, introduced the concept of the Fermi sea in a sem- 
inal 1926 paper lj. At low temperature all energy levels 



1 Historical note: Fermi surmised that these concepts would apply 
to bound electrons in an atom. The proper application to free 
electrons in a metal originated with Sommerfeld :2j. 
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single and double quantum dots has been studied exper- 
imentally [3 EH and theoretically [3 E E IM H3 ■ 
In a single quantum dot the entanglement can be be- 
tween the angular momentum of electron and hole form- 
ing a single exciton, but then the spatial separation of 
the electron-hole pair is problematic. In double quantum 
dots an exciton in one dot can become entangled with 
an exciton in the other dot, so the spatial separation is 
automatic. 

An altogether different theoretical proposal involves 
the "energy-time" entanglement of an electron in the con- 
duction band with a hole in the valence band, excited at 
a p-n junction by a laser with a long coherence time [22| . 
In this form of entanglement (originating from quantum 
optics the energy and time of creation of each par- 
ticle are uncertain, but the sum of the energies and the 
difference of the times arc well-determined. 



FIG. 2 Electron-hole pair (= exciton) formed by optical ex- 
citation of an electron from the (filled) valence band to the 
(empty) conduction band of a semiconductor. The key dif- 
ference with the electron-hole pair in the Fermi sea of Fig. 
is that the latter can be produced without optical excitation 
(because electron and hole coexist at the same energy, both 
within the conduction band). 

Kondo scattering by a magnetic impurity |l4j . 

Here we review the physics of entanglement in the 
Fermi sea from the three view points of production, de- 
tection, and utilization. The results presented in this 
review have all been published previously, except for the 
derivation in Sec. IIII.Bl of the critical temperature above 
which the electron-hole entanglement vanishes. 

In this introductory section we place our electron-hole 
entangler in a broader perspective by discussing similari- 
ties and differences with entanglers of excitons f Sec. ll.Jjf) 
and of photons (Sec. \l.C]i . 



B. Exciton entanglers 

In a semiconductor, holes can exist in the conduction 
band as well as in the valence band. A hole in the va- 
lence band bound to an electron in the conduction band is 
called an exciton, because it is produced by optical exci- 
tation of the electron across the band bap of the semicon- 
ductor (see Fig.|2J). In contrast, the electrons and holes in 
the conduction band of Fig. ^ coexist at the same energy, 
so that neither the creation nor the annihilation of the 
electron-hole pair involves any emission or absorption of 
a photon. The entanglement production and detection 
of electron-hole pairs in the Fermi sea is therefore purely 
electronic, without requiring any optical interface. 

Because the electron and hole that form an exciton 
are separated in energy by the band gap, their creation 
or annihilation is associated with photon emission or ab- 
sorption. Optically induced entanglement of excitons in 



C. Photon entanglers 

The production and detection of entanglement is much 
further developed in optics than in electronics. It is there- 
fore instructive to compare the mechanism for electron- 
hole entanglement of Fig. ^ with the two well-known 
mechanisms for photon entanglement illustrated in Fig. 
® 

The creation of a spin-entangled electron-hole pair by 
the applied voltage reminds one of the creation of a 
polarization-entangled photon pair by a pump laser in 
the process of spontaneous parametric down-conversion 
[24|. The key difference is that the electronic process 
is governed by a single-particle Hamiltonian, while the 
optical process relies on photon-photon interactions in a 
nonlinear crystal (in order to split a pump photon of fre- 
quency 2u> into two entangled photons of frequency lj). 
What the two processes do have in common is that the 
Bell pair is produced at random times, triggered by vac- 
uum fluctuations in the optical case and by stochastic 
tunnel events in the electronic case. The production rate 
is proportional to the nonlinear \^ coefficient and the 
tunnel probability r, respectively. Although the time of 
creation of a Bell pair is random, the two particles are 
precisely synchronized in both processes — without the 
need for any time resolution of the source. 

If we seek an analogy in the realm of linear optics, one 
can think of the entanglement of two photons of oppo- 
site polarization that are simultaneously scattered by a 
50/50 beam splitter [H |H HEl HJ . With probabil- 
ity 1/2 the two photons end up in different arms in a 
polarization-entangled Bell state. There are two require- 
ments for this process to succeed, which distinguish it 
from the electron-hole entangler: Firstly, noninteracting 
photons can not be entangled by a beam splitter if the 
sources are in local thermal equilibrium 29]. Somehow 
the Fermi sea, while being in local equilibrium, works 
around this optical no-go theorem. Secondly, the beam 
splitter relies on the indistinguishability of the two scat- 
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pump laser nonlinear crystal 



FIG. 3 Schematic illustration of the entanglement of electron-hole pairs by a tunnel barrier and of photons by a nonlinear 
crystal and a beam splitter. The differences and similarities of these entanglement mechanisms are discussed in the text. On the 
one hand, the Bell pair is produced spontaneously, without requiring synchronization at the source, by the voltage applied over 
the tunnel barrier and by the laser beam incident on the nonlinear crystal, while the beam splitter requires two synchronized 
single-photon sources. On the other hand, the tunnel barrier and the beam splitter both do not require interactions between 
the particles, while the nonlinear crystal requires photon-photon interactions to split the pump photon into two photons of 
lower frequency (down conversion). 



tered photons, which is only guaranteed if the incoming 
pair of photons are scattered simultaneously. (More pre- 
cisely, their wave packets should have a substantial over- 
lap at the beam splitter [3(J .) No such synchronization 
of the sources is required in the electronic case. 



II. ENTANGLEMENT BASICS 

The concept of entanglement has been made precise 
and quantitative in the context of quantum information 
processing |3l| . In this section we summarize what basic 
facts we need in order to calculate the amount of entan- 
glement produced by the application of a voltage over a 
tunnel barrier. 



A. Quantum versus classical correlations 

Loosely speaking, two spatially separated particles are 
entangled if their state can not be prepared from a prod- 
uct state by operating locally on each particle — not even 
with the exchange of classical bits of information. The 
elementary entangled state is the spin singlet 

|*bci 1 )-2- 1 / 2 (|T)aU) b -|I)a|T) b ), (2.1) 

also known as a Bell pair. The indices A, B label the two 
particles and the arrows f, J. indicate two spin states. The 
state I^Bcii) can not be prepared locally starting from a 
product state such as |T)a|4)b- 

The spins in the Bell state are correlated; a measure- 
ment of one spin, say with outcome |, projects the other 
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spin on the opposite spin state J.. The same applies to 
the mixed state with density matrix 

ft-ixed = i(iTXTi) A ®(u)ai) B +i(u>ai)A®(it>ai) B . 

(2-2) 

The difference between the quantum correlation of the 
pure entangled state 1|2.1|) and the classical correlation of 
the mixed nonentangled state l|2.2[) is that the quantum 
correlation persists if we measure the spin along a dif- 
ferent axis, while the classical correlation is diminished. 
For example, the change of basis 

|T)-2- 1 / 2 (|t) + U)), |l)-2- 1 / 2 (|t)-||)) (2.3) 

(both for spin A and spin B) leaves I^Beii) invariant, but 
Pmixed becomes a mixture of parallel and anti-parallel 
spins. 

B. Bell inequality 

The Bell inequality is a test to distinguish quantum 
from classical correlations by comparing the correlation 
along different directions ■ If spin A is measured along 
unit vector a and spin B along unit vector b, then the 
correlator is the expectation value 

C ab = {(a-(r) A ®(b-<T) B ), (2.4) 

with er = (cr x) u y ,(J z ) the vector of Pauli matrices 

^ = (J J) . *v = (? V) ' *' = (o -l) ' (2 ' 5) 

The Clauser-Horne-Shimony-Holt (CHSH) form of the 
Bell inequality reads [33| 

B=\C ab + C a , b + C aV -C a , v \<2, (2.6) 

for classically correlated spins. If B > 2 for some choice 
of unit vectors a,b,a'b', then the spins are quantum 
correlated. The converse of this statement is not true: 
There exist mixed states that are entangled (in the sense 
that they can not be prepared locally) yet satisfy B < 2 
for all sets of unit vectors (see Sec. III.Dll . 

An advantage of the CHSH form of the Bell inequality 
is that it can be used even if the detection efficiency is 
smaller than 1. If a coincidence detection of spins A and 
B finds N± G {0, 1} spins pointing up or down at each 
detector, then the Bell-CHSH inequality reads 

B = \((N+(a)-N-(a))(N + (b)-N-(b))) 

+ ({N + (a') - N-(a'))(N+(b) ~ N-(b))) 
+ ((N + (a) N.(a))(N + (b') - iV_ (&'))) 

({N + (a') - N_(a')){N + (b') - N.(b')))\ < 2u, 

(2.7) 

v = ((N+{a)+N-(a))(N+(p)+N-(b))), (2.8) 

where the detection efficiency v is assumed to be inde- 
pendent of the unit vectors a and b along which spins A 
and B are measured. 



C. Entanglement measures for pure states 

It is rather straightforward to quantify the amount of 
entanglement present in a pure state, so we discuss this 
case first — before proceeding to entanglement measures 
for mixed states in the next subsection. 

Since an entangled state can not be prepared locally, 
one needs to exchange a certain amount of quantum infor- 
mation to create it out of a product state. This quantum 
information can take the form of Bell pairs, shared be- 
tween A and B. Bell pairs play the role of a "currency" 
(1 Bell pair = 1 bit of entanglement = 1 ebit), by means 
of which one can quantify entanglement. The average 
number of Bell pairs per copy needed to prepare a large 
number of copies of the pure state |\&) is given by 

£ = -Ti A p A \og 2 p A , p A =Tr B \*)(*\. (2.9) 

(The reduced density matrix p A of subsystem A is ob- 
tained by tracing out the degrees of freedom of subsystem 

The quantity £ is called the entanglement entropy, or 
entanglement of formation. It can be used to quantify the 
amount of entanglement present in the pure state of any 
system with a Hilbert space H = H A ® TLb that can be 
divided into two subsystems A and B (a socalled bipartite 
system). It is important to emphasize that entanglement 
entropy is defined relative to a particular partitioning. 
Indeed, a system might well be separable relative to Ti = 
Tt A ®H.b and entangled relative to H = Ti.' A ® Ti! B . 

In the special case that each subsystem contains a sin- 
gle qubit, one can use an equivalent measure of entan- 
glement called the concurrence. A normalized two-qubit 
state |\&) has the form 

2 2 

l*)=EE^I^W fl ' Tr 77 t = l, (2.10) 

where we have labeled 1 =|, 2 =|. Substitution into Eq. 
(|2.9|) gives the following expression for the entanglement 
entropy in terms of the 2x2 matrix of coefficients 7: 

£ = J r (i + i v / l-4Dct77t), (2.11) 

where the function T(x) is defined by 

T{x) = -x\o% 2 x- (l-x)log 2 (l-a;). (2.12) 
The quantity 

C = 2-v/Det77t (2.13) 

is called the concurrence of the two qubits. It is in one- 
to-one relation with £ but has a somewhat simpler ex- 
pression. A Bell pair has unit concurrence and carries 
one bit of entanglement, while both C and £ vanish for a 
product state. 

From an experimental point of view, the Bell inequality 
is the most direct way to quantify the degree of entan- 
glement of two qubits. By maximizing the correlator B 
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in Eq. I|2.6[l over the unit vectors one measures the Bell 
parameter 

S max = max B. (2.14) 

a.b,a' b' 

The relation between B maK and C is 34] 

S max - 2Vl + C 2 . (2.15) 

The Bell parameter varies from 2 to 2^2 when the con- 
currence varies from to 1 . As mentioned in the previous 
subsection, any entangled pure state violates the Bell in- 
equality ($ max > 2 if and only if C > 0) . 

D. Entanglement measures for mixed states 

The density matrix p of a mixed state can be decom- 
posed into pure states \^ n ) with positive weight p ni 

n n 

The states l^n) are normalized to unity (^l^n) = 1, 
but they need not be orthogonal. The convex-sum de- 
composition i|2.16[) is therefore not unique — there are 
many equivalent representations of p as a mixture of pure 
states. 

A mixed state in the bipartite Hilbert space Ha <8> Ti-B 
is nonentangled (= separable) if there exists a convex- 
sum decomposition <|2.16[) into pure product states |35j . 
meaning that |* n ) = |$„)a|$^)b with \$ n ) A € U A and 
|$^)b € Hb for all n. The entanglement entropy £ of a 
separable mixed state vanishes, with the definition 

£ = min VV£(*„). (2.17) 

{*»,p»} v 

Here £ (^n) is the entanglement entropy of the pure state 
l^n), defined in Eq. (|2.9p . and the minimum is taken over 
all convex-sum decompositions <|2.16ll of p. 

The entanglement entropy, or entanglement of forma- 
tion, of a mixed state is an upper bound to the average 
number of Bell pairs per copy that it costs to prepare 
many copies p®p®p<®- ■ ■ of the state p j3|J. This socalled 
entanglement cost would be equal to £ if £ would be an 
additive quantity meaning that £{p® p) — 2£{p). The 
additivity of the entanglement entropy has been proven 
for certain classes of mixed states [37| , but not in general. 
It might be that for certain mixed states entanglement 
has a "discount" , in the sense that the cost per copy is 
less for many copies than for a single copy. 

The entanglement entropy is also an upper bound to 
the average number of Bell pairs per copy that one can 
extract (or "distill") from many copies of an entangled 
state, using only local operations and classical commu- 
nication |38j |. The production of entangled pure states 
out of Bell pairs is reversible, meaning that £ equals the 
distillable entanglement. For mixed states this is not the 



case in general: The distillable entanglement can be less 
than £, so a fraction of the Bell pairs can be lost in the 
conversion from Bell pairs to entangled mixed states and 
back to Bell pairs. 

The calculation of £ for a general multi-qubit bipartite 
mixed state is more complicated than for a pure state, 
because of the need to minimize over convex-sum decom- 
positions. In the case of two qubits a closed form expres- 
sion for £ exists, due to Wootters [39| : 

£ = ^(i + i V / l-C 2 ), (2.18) 
C = max{0, v/aT - ^/M - \/% - y/M}. (2.19) 
The Xi 's are the eigenvalues of the matrix product 

P- (<Jy<g) Gy) ■ P* ■ ((Ty ® <7y), 

in the order Ai > A2 > A3 > A4. The quantity C is again 
called the concurrence. For a pure state (with p 2 = p) 
the definition (|2.19() of C is equivalent to Eq. I|2.13|l . 

As a simple example we take the rotationally invariant 
mixed state 

|0Werner(O - J (1 - t)U + £|*Bell) (*Bcll I , ~| < f < 1, 

(2.20) 

with I4 the 4x4 unit matrix. This is the socalled Werner 
state |35j. The concurrence is 

C = max{0, i(3£-l)}. (2.21) 

The Werner state is separable for |£| < 1/3. Notice that 
separable mixed states, unlike separable pure states, oc- 
cupy a region in parameter space of nonzero measure. 

For a mixed state there is no one-to-one relation be- 
tween the Bell parameter (|2.14|) . which quantifies the 
maximal violation of the CHSH inequality l|2.6|l . on the 
one hand, and the entanglement entropy or concurrence, 
on the other hand. Depending on the density matrix, 
B max can take on values between 2C\/2 and 2Vl + C 2 . 
The general formula 0j [13 

# max - 2y/ Ul + u 2 (2.22) 

for the dependence of S max on p involves the two largest 
eigenvalues ui, of the real symmetric 3x3 matrix R T R 
constructed from Rki — Tr p Ok ® <Ji ■ In the case of the 
Werner state (I2.2U|) one finds £> max = 2£y / 2, so this state 
is entangled without violating the CHSH inequality if 
1/3 < £ < 1/V2. 

E. Particle conservation 

As discussed in Sec. III.AI to distinguish classical from 
quantum mechanical correlations it is essential that the 
strength of the correlation is compared in different bases. 
In other words, to detect entanglement between two sep- 
arated qubits one needs to be able to rotate each qubit 
individually creating superpositions like Eq. I|2.3|) . As 
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emphasized by Wiseman and Vaccaro |42j, if a conser- 
vation law prevents the rotation, then the entanglement 
reduces effectively to a classical correlation — becoming 
inaccessible as a quantum resource. 

For electrons the restriction on the accessible entan- 
glement imposed by conservation of particles is of pri- 
mary importance. Particle conservation prevents the lo- 
cal creation of a state that is in a superposition of a 
different number of particles: If |0) denotes an empty 
single-electron state and |1) denotes a filled state, then 
the superposition 

|0)~ 2-^(10) + |l)), 2-^(10) (2.23) 

can not be created locally. As a consequence, the single- 
electron state 

|* 1 )=2- 1 /2(| ) A |1) S -|1) A |0) B ) (2.24) 

has £ = 1, like the two-particle Bell state 1)2. but this 
entanglement is inaccessible because of particle conserva- 
tion. (Such useless entanglement has been dubbed "fluffy 
bunny" entanglement 43].) 

To account for particle conservation in a bipartite 
multi-electron state we need to project the density ma- 
trix p onto sectors N pq of Fock space with an integer 
particle number p at A and q at B. This projection elim- 
inates unobservable coherences between sectors of differ- 
ent particle number. The total accessible entanglement 
entropy £ pa rt, constrained by particle conservation, is the 
weighted sum of the entanglement from the individual 
sectors [42] |. 

£ P art = y^,w pq £(p pq ), (2.25) 

Ppq = RpqPttpq, Wpq = ^^pqP^pq- (2-26) 

Here TL pq is the projection operator onto N pq . 

For the state (|2.24|) we have two sectors p = 0, q = 1 
and p = l,q = 0, with wqi — wiq — 1/2 and poi — 
(|0)(0|)4|1>(1|) B) Pio = (|1)(1|)a(|0)(0|) b . The accessi- 
ble entanglement 

£ P art = \£{pm) + \£{pva) = ( 2 - 27 ) 

vanishes, as expected, because poi and pio are product 
states. 



F. Phase reference 

The restriction on the accessible entanglement imposed 
by particle conservation can be (partially) removed if the 
two parties A and B share a phase reference |44j, |45j, |4^| . 
In the quantum optical context, a phase reference is an 
identical pair of coherent states (one at A and one at B). 
In the electronic context of interest here, the most natural 
candidate for a phase reference is a pair of superconduc- 
tors with the same phase of the order parameter (which 



we may therefore set equal to zero). Andreev reflection 
at the superconductor converts a state with N electrons 
into a state with N ± 2 electrons, adding or extracting a 
Cooper pair to or from the superconducting condensate. 
By combining Andreev reflection and normal reflection 
one can realize the superpositions 

|0>^2- 1 / 2 (|0) + |2)), |2)-2- 1 / 2 (|0)-|2)). (2.28) 

This is indeed how the charge qubit in a socalled "Cooper 
pair box" is rotated |47| . 

The rotation (|2.23|) . involving a charge difference of le 
rather than 2e, remains inaccessible, in view of the lack 
of a macroscopically coherent state of unpaired electrons. 

III. HOW TO ENTANGLE FREE PARTICLES 

One way to entangle particles is by letting them in- 
teract with each other. If the Hamiltonian H = Ha + 
Hb + Hab of two spatially separated particles A and B 
contains an interaction term Bab-, then their time de- 
pendent state e~ lHt l h Y5>) aV&) b will evolve from a prod- 
uct state into an entangled state. Since the particles are 
spatially separated, and hence distinguishable, it is irrel- 
evant whether they are fermions or bosons. 

Indistinguishable particles may become entangled by 
scattering from an external potential, such as the beam 
splitter or tunnel barrier of Fig. even if they do not 
interact with each other. A fundamental difference now 
appears between fermions and bosons, as we discuss in 
the following subsections. 

A. Free bosons 

Whether or not a beam splitter is able to entangle 
noninteracting particles depends on the quantum state 
of the sources [2^, l26l 1271 12H 129^ . Free bosons can not be 
entangled if the sources are in (local) thermal equilibrium. 
For a proof of this "no-go theorem" we follow Xiang-bin 

We consider the geometry of Fig. 0] The density ma- 
trix p m of the multi-mode state incident from the sources 
onto the beam splitter can be written in the coherent 
state representation, 

p hl = J daP(a)\a){a\, (3.1) 

\a) =exp(a t -a-a-a*)|0). (3.2) 

The coherent state |a) = |ai, 02, . . .) is an eigenstate of 
the annihilation operator a n with complex eigenvalue a n . 
We have abbreviated da = FT <£Rc a n dlma n . 

The real function P(a) may take on negative values. 
In that case the state p m is called "nonclassical" , because 
P{ot) can not be interpreted as a classical probability 
distribution [24|. A thermal state is a classical state. (It 
has a Gaussian P(a) oc exp(— J2 n |a„| 2 //„), with /„ the 
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FIG. 4 Schematic setup for entanglement of indistinguishable 
particles by a beam splitter. Two sources Sa and Sb create 
a separable initial state. The beam splitter scatters the parti- 
cles into two detectors Da and Db- An entangled final state 
may result if the detectors can not distinguish from which 
source a particle originated. For bosons no entanglement can 
be produced by the beam splitter if the sources are in ther- 
mal equilibrium. For fermions this "no-go theorem" does not 
apply. 



Bose-Einstein occupation number of mode n.) Nonclas- 
sical states include photon-number (or Fock) states and 
squeezed states. 

The beam splitter transforms the annihilation opera- 
tors a — {ai, et2, . . .} of the incoming state into annihi- 
lation operators b = {t>i, &2, - ■ ■} of the outgoing state. 
This is a unitary transformation, 



S 



nm^m j 



s- s f 



l, 



(3.3) 



defined by the scattering matrix S. The density matrix 
Pout of the outgoing state is obtained from p m with the 
help of Eq. {3~3Jl. The result is 



Pout= / d/3P(S^(3)\/3){l3\, 



|/3)=cxp(bt./3-6./3*)|0). 



(3.4) 
(3.5) 



We have made a change of variables from a to (3 = S ■ a 
and used that da. = d(3 because of the unitarity of S. 

We now group the outgoing modes into Na modes that 
are detected by A and Nb modes that are detected by B. 
Since the creation and annihilation operators of different 
modes commute, we may write \(3) as a product state, 



1/3} 



/ N A 

n 

\n=l 



/N A +N B 

n 

\n=N A + l 



bl/3„-b n l3* 



|0>. 



(3.6) 

The pure state \(3) is therefore separable for all (3. Refer- 
ring to the definition of a separable mixed state in Sec. 
III.DI we conclude that p out in Eq. (|3.4I) is separable if 
P > for all (3, because in that case the density ma- 
trix has a decomposition into separable pure states with 
positive weights. 



Nonclassicality of the sources is therefore a necessary 
condition for entanglement by single-particle scattering 
[2fj| . This condition is violated by a thermal source, prov- 
ing the no-go theorem. 



B. Free fermions 

In contrast to bosons, fermions can be entangled by 
single-particle scattering even if the sources are in (lo- 
cal) thermal equilibrium The no-go theorem of the 
previous subsection does not apply, because it relies on 
a distinction between classical and nonclassical bosonic 
density matrices without a fermionic analogue. 2 To 
demonstrate the free-fermion entanglement in its sim- 
plest form, we consider a single-mode, 3 spin-degenerate 
conductor at zero temperature, containing a barrier with 
spin-independent tunnel probability r. 

We refer to the geometry of Fig. ^ Tunnel attempts 
in the energy range eV occur with a rate eV/h. indepen- 
dently for spin- up and spin-down electrons |49j| . In each 
time interval h/eV, either 0, 1, or 2 electron- hole pairs 
are created, transforming the unperturbed Fermi sea |0) 
into the superposition 4 



10) 



(l-T)\ 0)-e 2 ^T \n)H\U)e 

- e i V2r(l-r)2- 1 /2 (| T ) fc | T)e + | l)h \ i)e ), 



Bell pair 



(3.7) 



with <j) the phase difference of the (spin-independent) 
transmission and reflection amplitudes. In each prod- 
uct of kets |-), the first ket refers to the hole excitations 
at the left of the barrier and the second ket refers to the 
electron excitations at the right of the barrier. The third 
term (labeled "Bell pair" ) involves the creation of a sin- 
gle electron- hole pair (with spin up or spin down), while 
the first and second term involve the creation of and 2 
electron- hole pairs, respectively. 

As discussed in Sec. lILEl the total entanglement £ par t, 
constrained by particle conservation, is the sum of the en- 
tanglement carried by the 0, 1, and 2-particle states. This 
gives (1 - t) 2 x + 2r(l - r) x 1 + t 2 x = 2t(1 - t) bits. 



2 The fermionic analogue of the coherent state decomposition 13.11 
involves an integr atio n over Grassman variables, rather than over 
complex numbers |4S| . A Grassman integral is not a valid convex- 
sum decomposition of a density matrix, because it lacks a positive 
integration measure. 

3 The single-mode restriction is not essential in the case of spin- 
independent scattering; in the eigenbasis of the transmission ma- 
trix product ttt, multiple modes represent multiple "copies" of 
the spin-entangled state, so the total entanglement production is 
simply the sum of the contributions from the individual modes. 

4 Tunneling creates a coherent superposition rather than an in- 
coherent mixture because of the indistinguishability of the elec- 
trons. 
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To obtain the entanglement production in a detection 
time tdet we should multiply by the number i<j ct x eV/ h 
of tunnel attempts in that time interval, resulting in 

2eV 

£ P artAdot = —r— r(l - r). (3.8) 

The maximum entanglement production rate of eV/2h 
is reached for r = 1/2 (corresponding to a 50/50 
beam splitter). In the tunneling limit r <C 1 one has 
£ P art/*det = 2eVr/h. 

Notice that the states at the left-hand-side and at 
the right-hand-side of Eq. I|3.7I) are related by a single- 
particle unitary transformation, and yet the left-hand- 
side is separable while the right-hand-side is entangled. 
There is no contradiction because entanglement is only 
invariant under local unitary transformations. Scattering 
by the tunnel barrier mixes the degrees of freedom at the 
left and at the right, so it is a nonlocal transformation — 
which has the capability to transform a separable state 
into an entangled state. 



The result (|3.8(l points to an important connec- 
tion between entanglement production and shot noise: 
The entanglement production rate in the case of spin- 
independent scattering is proportional to the spectral 
density P no ise of the current fluctuations through the con- 
ductor, given by (H3> E3 

2e 2 

P noisc = 2eV—r(l -t) = 2e 2 £ part /i dot . (3.9) 

We will return to this connection in Sec. IVIII 

In Ref. the more general case of spin-dependent 
scattering was analyzed (still for a single-mode conduc- 
tor and at zero temperature) . The calculation is reviewed 
in Appendix 1X1 The transmission matrix t of the tunnel 
barrier is a 2 x 2 complex matrix. The transmission eigen- 
values t\,T2 € [0, 1] are the eigenvalues of the Hermitian 
matrix product tt' . The result for the entanglement pro- 
duction rate is 



-'part 



/tde 



eV 
~~h 

eV 



(n + t 2 - 2t x t 2 )T (± + yi-C 2 



{(ti + t 2 - 2tit 2 ) log 2 (ri + r 2 - 2tit 2 ) 



- n(l - r 2 )log 2 [ri(l - r 2 )] - r 2 (l - n) log 2 [r 2 (l - n)]} 
2Vrir 2 (l-Ti)(l-r 2 ) 
n + r 2 - 2tit 2 



(3.10) 
(3.11) 



The concurrence ()3.11(l of the electron-hole pair reaches 
its maximum value C = 1 for the spin-independent case 
T\ = t 2 = r. In that case Eq. (|3.10(l reduces to Eq. JIOJ. 

Note that the concurrence is entirely determined by 
the transmission eigenvalues; the eigenvectors of tv do 
not contribute. This means, in particular, that chan- 
nel mixing does not degrade the entanglement as long 
as the transmission eigenvalues remain unaffected. Note 
also that the direct relation (|3.9|l between entanglement 
production rate and shot noise power no longer exists 
in the case of spin-dependent scattering. The shot noise 
power in this case equals P no ise = 2eV (e 2 / '/i)[ri(l — n) + 
t 2 (1 — r 2 )], which is not in one-to-one relationship with 
Eq. {EBB . 

So far we have assumed zero temperature. At nonzero 
temperature the electronic state is mixed rather than 
pure, and the rate of entanglement production decreases. 
General considerations require that the entanglement 
vanishes identically for temperatures T greater than a 
critical temperature T c . The calculation of T c is pre- 
sented in Appendix [51 for the case of spin-independent 
scattering. The result (|B18|) is plotted in Fig. El We find 



that T c is given by 



r(l - r)sinh 2 (eF/2fc s T c ) = 1/4, 



(3.12) 



as plotted in Fig. EI When T c is approached from below, 
the entanglement vanishes with a logarithmic singularity, 

£ P art oc (T c - T f ln(T c - T)-\ when T T T c . (3.13) 

The maximal critical temperature for a given voltage is 
reached for r = 1/2, when T c = 0.57 eV/ks- In the 
tunneling regime one has 5 



eV 



k B ln(l/r) 



, if r < 1. 



(3.14) 



Because the dependence of T c on r is only logarithmic, 
the critical temperature will be of order eV/ks even for 
very small tunnel probabilities. 



A second critical temperature T' c exists, above which the Bell- 
CHSH inequality can not be violated. In general T' c < T c . In the 
tunneling regime r<l the two critical temperatures coincide. 
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FIG. 5 Temperature dependence of the entanglement pro- 
duction £ P art (in units of So = eVtdet/h), calculated from 
Eq. (IB18pi . Results are shown for three values of the spin- 
independent transmission probability r. (The same plots are 
obtained if r is replaced by 1 — r.) 
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FIG. 6 The entanglement production of Fig. is identically 
zero for T > T c , with the critical temperature given by Eq. 
13.12H and plotted in this figure. 



We close this subsection by calculating how much the 
constraint of particle number conservation reduces the 
entanglement production £ par t below the unconstrained 
value £ . We take zero temperature, so we can use the 




FIG. 7 Comparison of the unconstrained entanglement pro- 
duction S [given by Eq. 13.15H 1 with the entanglement pro- 
duction fpart constrained by particle conservation [given by 
Eq. 13.101 1. Both quantities are normalized by So = eVtdet/h 
and plotted as a function of the spin-independent transmis- 
sion probability ri = T2 = r. Zero temperature is assumed. 



pure-state formula i|2.9|) for £. The result is 6 

£/t d et = ^[Hn)+HT2)], (3.15) 

which in the tunnel limit simplifies to 

eV 

S/tdet = -r{-n log 2 Ti - r 2 logr 2 ), if Ti, t 2 < 1. 
h 

(3-16) 

The constrained entanglement production (|3.1Ufl simpli- 
fies in the tunnel limit to 

eV 

£part Adet = — [~Tl log 2 T\ - T 2 log T 2 

+ (ti + r 2 ) log 2 (ri + n)], if n,T 2 < 1. 

(3-17) 

We conclude that the constraint of particle conservation 
reduces the entanglement production rate by an amount 
(eV/h)(ri +r 2 ) log 2 (ri +r 2 ) in the tunnel limit. The two 
quantities £ and £ par t are compared in Fig.[7|for the case 
of spin-independent transmission probability. 



IV. SPIN VERSUS ORBITAL ENTANGLEMENT 

Any normal-metal conductor containing a localized 
scatterer can be used to entangle the outgoing states to 
the left and right of the scatterer. The particular imple- 
mentation described originally in Ref. [3j uses tunneling 
between edge channels in the integer quantum Hall ef- 
fect. A variety of other implementations have since been 
proposed, as shown in Figs. 0D and [5] 

In all implementations a stream of entangled electron- 
hole pairs is produced by application of a voltage over 
a normal-metal conductor. What differs is the degree 
of freedom which is entangled and the means by which 
the entanglement is detected. The three implementa- 
tions from Refs. jH^, L3. EH shown in Fig. entangle 
the spin degree of freedom and use ferromagnetic con- 
tacts to measure the spin-resolved correlator needed to 
test for violation of the Bell inequality (as proposed by 
Kawabata j5|| and Chtchelkatchev et al. |57|). The three 
implementations from Refs. |3. l58ll59]| shown in Fig. [^en- 
tangle a spatial (= orbital) degree of freedom instead of 
spin and use non-ferromagnetic contacts for the detection 
jas^ proposed by Samuelsson, Sukhorukov, and Biittiker 

The implementations of Fig.[S]have the advantage that 
a spin-entangled electron-hole pair is much less sensitive 
to decoherence than an orbitally entangled pair, simply 
because most external degrees of freedom couple to the 



6 The entanglement production rate given without further discus- 
sion in footnote 24 of Ref. |3| is this unconstrained value £, rather 
than the more appropriate constrained value £ pa rt- 
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FIG. 8 Three geometries to produce and detect spin entanglement in a normal conductor with ferromagnetic contacts. Panels 
a, b, and c are taken, respectively, from Refs. |53j . [54) . and [5Sj|. The magnetization axis in the contacts is indicated by the 
vectors a, b or by arrows. 



position of the electron rather than to its spin. The en- 
tanglement production is optimal because, assuming that 
the scattering by the normal-metal conductor is spin- 
independent, each electron- hole pair represents a max- 
imally entangled Bell pair (= spin singlet). The maxi- 
mum entanglement production rate of eV/2h bits/second 
is reached for a 50/50 beam splitter [cf. Eq. (|3.8|) ]. There 
is no need for the conductor to be single-mode — multi- 
ple orbital modes simply produce multiple copies of the 
spin-entangled Bell pair. The disadvantage of these im- 
plementations is that it is difficult to resolve the different 
spin components of the current. A ferromagnetic contact 
with an adjustable magnetization axis may provide the 
required spin filtering capability 60]. In order to avoid 
backscattering, the conductance of the non-ferromagnetic 
conductor should be much less than the total conduc- 
tance of the ferromagnetic contacts [55|. An alterna- 
tive non-ferromagnetic spin filter might be provided by a 
quantum dot in a strong magnetic field |6ll |62j . 

The orbital entanglement in Figs. [5Kb, c refers to differ- 
ent spatial degrees of freedom. The quantum dot of Fig. 

has two single- channel point contacts at the left (la- 
beled Li, L 2 ) and two more at the right (labeled i? 2 ). 
Each point contact plays the role of one of the spin direc- 
tions in the spin entanglers. In the Hall bar geometry of 
Fig. |5J} this role is played by the two quantum Hall edge 



channels (labeled by their Landau level index 1,2) at the 
left and right of the constriction. There is only a single 
Landau level in the disk geometry of Fig. [U — the role 
of the spin direction being played by the source contact 
2 or 3 from which an electron reaches the detectors at A 
and B. 

The quantum dot entangler of Fig. [§]} has the advan- 
tage that no magnetic field is required for its operation, 
but the disadvantage that backscattering at the detector 
(back into the quantum dot) can not easily be avoided. A 
strong perpendicular magnetic field, as in Figs.[5|D,c, sup- 
presses backscattering: The Lorentz force constrains the 
electrons to move unidirectionally along equipotentials at 
the edge of the two-dimensional electron gas. (For a tu- 
torial on edge channel transport, see Ref. [63 •) Let us 
discuss these two edge channel entanglers in some more 
detail. 

The Hall bar entangler contains a beam splitter formed 
by a split gate electrode (shaded rectangles at the center 
of Fig. I^Jd), as realized experimentally in Ref. .64]. The 
current in edge channel 1 can be measured independently 
of that in edge channel 2 by diverting an equipotential 
into a point contact, as in Ref. 65]. This capability to 
filter edge channels is not enough for a test of the Bell 
inequality; the filter should also have an adjustable eigen- 
basis. In the spin filters of Fig.^Jthe eigenbasis is just the 
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FIG. 9 Three geometries to produce and detect orbital entanglement in zero magnetic field (a) and in the quantum Hall effect 
regime (b,c). Panels a, b, and c are taken, respectively, from Refs. |H^|. |!J, and |H?J. The solid (dashed) lines in panels b,c 
indicate filled (empty) edge channels, with the arrows pointing in the direction of motion. 
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magnetization axis, which can be rotated by a magnetic 
field. Equivalently, one could rotate the spin at fixed 
magnetization axis of the spin filter. The edge channel 
filters have a fixed eigenbasis, so one would need to "ro- 
tate" or mix the edge channels before they are filtered. 
The two gate electrodes labeled Ul and Ur in Fig. Et> 
should serve this purpose, by creating a strong local elec- 
tric field that causes inter-edge-channel transitions. Such 
edge channel mixers have not yet been demonstrated in 
the laboratory, however. 

The electron-hole pairs produced in this geometry have 
a concurrence given by Eq. (|3.11|) . For an optimal en- 
tanglement the two transmission eigenvalues t\ and Ti 
should not be widely different. This is problematic, since 
edge channel number 1 tunnels over a shorter distance 
through the constriction than edge channel number 2, 
and therefore T2 will in general be <C n . It may be pos- 
sible to resolve this problem by introducing disorder into 
the constriction, in order to mix the edge channels while 
they are tunneling and equalize the transmission eigen- 
values. 

Samuelsson, Sukhorukov, and Biittiker [59] found a 
way to avoid the need to mix edge channels from different 
Landau levels by working with a single Landau level in a 
disk geometry (Fig. . This geometry is the electronic 
analogue of the Hanbury-Brown-Twiss interferometer in 
quantum optics |24| . Edge channels circulate counter- 
clockwise along the outer edge of the disc and clockwise 
along the inner edge. Four split-gate electrodes (labeled 
A,B,C,D) couple the inner and outer edge. Four electrical 
contacts (numbered 2,3,5,8) are connected to the outer 
edge and four more (numbered 1,4,6,7) are connected to 
the inner edge. All contacts are grounded, except con- 
tacts 2 and 3 at voltage V. This number 2,3 of the source 
contact replaces the Landau level index 1,2 in the Hall 
bar geometry. The split gates A and B replace the edge 
channel mixers Ul and Ur, while the two split gates C 
and D together replace the beam splitter at the center of 
the Hall bar. By adjusting the gates C and D separately 
one can readily equalize the two transmission eigenvalues 
and ensure that the electron- hole pair arriving at gates A 
and B is maximally entangled. Tunneling through split 
gates A and B provides the mixing needed to test the 
Bell inequality Such mixers within a single Landau level 
have been realized experimentally j6|J . 



V. ENTANGLEMENT DETECTION BY NOISE 
MEASUREMENTS 

A. Tunneling regime 

We emphasized in our presentation of the electron- 
hole entangler, in Sec. ^ that no synchronization at the 
sources is required: The electron and hole are automat- 
ically produced at the same instant in time, since they 
originate from a single tunnel event. We did not yet 
address the issue whether synchronization might be re- 
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FIG. 10 Current pulses of electrons (solid) and holes (dashed) 
produced at a tunnel barrier with transmission probability 
r« 1. The width St of the pulses is smaller than their spacing 
At by a factor r. Low- frequency noise measurements correlate 
the spin of each electron with that of every hole, but only the 
spins from coincident current pulses contribute (because the 
spins from non-coincident pulses are uncorrelated) . In this 
way the electron-hole entanglement can be detected without 
requiring time-resolved measurements. 



quired at the detectors. 

In a typical optical experiment [24) . the two photon 
counters that test for entanglement of a photon pair are 
synchronized by a coincidence circuit. The time reso- 
lution required to detect the current pulses from indi- 
vidual tunnel events is difficult to reach in solid state 
electronics. Fortunately, it is not needed. As pointed 
out by Samuelsson, Sukhorukov, and Biittiker [I3j . low- 
frequency noise measurements can test for the violation 
of the Bell inequality even if individual current pulses are 
not resolved. 

A sequence of current pulses through a tunnel barrier 
is shown schematically in Fig. ^| Their characteristic 
width St = h/eV is much less than their mean spacing 
At = e/I = h/2erV if the tunnel probability r < 1. It 
is therefore possible in principle to detect the individual 
entangled electron-hole pairs and test for a violation of 
the Bell inequality by correlating their spins. In a low- 
frequency noise measurement a large number of current 
pulses is detected and the spin of each electron is corre- 
lated with the spin of every hole, regardless of whether 
the two particles were produced by the same tunnel event 
or not. Electrons and holes from different tunnel events 
have uncorrelated spins, so they do not contribute to the 
correlator. The low-frequency measurement is therefore 
equivalent to a time-resolved measurement. 

As formulated in Sec. III.BI the CHSH inequality con- 
tains the two-particle correlator C a b of spin A mea- 
sured along unit vector a and spin B measured along 
unit vector b. In a conductor one would measure cor- 
relators of currents rather than of individual particles. 
Chtchelkatchev, Blatter, Lesovik, and Martin [57j have 
reformulated the Bell inequality in these terms. 

The correlator Kij = (Na^Nbj) of the number of elec- 
trons (holes) detected in a time t^ct with spin i (j) at 
detector A (B) is obtained by integrating the current 
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correlator 
Ku = 



tdet pt-det 

dt dt' (I Al (t)I BJ (t')) 

n Jo 

2 



tdet^iX/Bj) + 



WMdet/2) 
^ U ' ( )_ 

(5.1) 



Here Cjj(w) is the frequency dependent correlator of cur- 
rent fluctuations I67I. 



This non-unit detection efficiency appears in the denom- 
inator of Eq. (|5.3I) , so that the resulting correlator counts 
only the successful measurements [in accordance with Eq. 

one hand, tdct St = 

large one may approximate 
idct^(^) in Eq. I|5.1|) by its long- 



Because, on the 
eV/h is sufficiently 
2 sin 2 (wtdet /2)/ttuj 2 
time limit, hence 



A", 



t deti 1 A,i) (I B,j) +idetCy(0). 



(5.4) 



, l-|^/eF| if |7iu;/eF|<l, 
■ 1 1 if |/iw/eV| > 1. 



(5.2) 

We have assumed that V is bigger than temperature 
but still small enough that the energy dependence of 
the scattering matrix may be neglected in the range 
(Ep, Ep + eV). The dependence of Kij on the unit vec- 
tors a, b (along which the spin is measured) is implicit. 
The subscripts j,je{t,j}={l,2} may equally well re- 
fer to two orbital degrees of freedom, as explained in Sec. 

mi 

The correlator C a b in the CHSH inequality Q2.6|l is 
determined by 



Cab — 



((N Ai1 -N Ai2 )(N B s-N b ^)) 
{(N AA +N At2 )(N BA + N B , 2 )} 
if n + K 22 - K 12 - K 21 

Kn + K 22 + if 12 + K 2 i ' 



(5.3) 



where the detection time tdct should be in the range 
St <C tdot <^ At to ensure coincidence detection. The 
average (•••) is over many measurements, of which a 
fraction v = ((N AA + N A ^ 2 )(N BA + N B ^ 2 )) is success- 
ful in detecting one particle at each detector A and B. 



J 



On the other hand, because tdct <C At = h/2erV is 
sufficiently small one may neglect the first term on the 
right-hand-side of Eq. (|5.4|l relative to the second term: 
tdct(I) 2 /C ~ tdet(-0/e <C 1. The correlator C a {, in the 
CHSH inequality is then given entirely in terms of low- 
frequency current correlators, 



C a 



Cn(0) + C 22 (0)-C 12 (0)-C 21 (0) 



(5.5) 



Cu (o) + c 22 (o) + C12 (o) + c 2 i (o) • 

To verify this line of argument, one can calculate the 
concurrence of the electron-hole pairs that follows from 
the violation of the CHSH inequality using the correlator 
(|5.5() and compare with the expected result (|3.11() . This 
calculation [3j is outlined in Appendix lC.il 



B. Beyond the tunneling regime 

If the transmission probability t is not <C 1 , then there 
is no separation of the time scales St, At and the argu- 
ment from the previous sub-section breaks down. It is 
still possible to determine the degree of entanglement us- 
ing low-frequency measurements, but a modified expres- 
sion for the spin-spin correlator C a b is needed [5^, . 

For short detection times tdct <C h/eV, shorter than 
both St and At, one may take the limit 



2sin>t d et/2) ^ 
dui 7, Uij{oj) = 



duj 



eV 



27r -a» = — Cy(0), 



cf. Eq. (|5.2|) . Instead of Eq. (|5.4|l we now have 

Kij - t 2 dct [(I Al )(I Bj ) + (eV/KjdM]- 
Substitution into Eq. I|5.3|l gives the spin-spin correlator 

{h/eV)(I AA - I A , 2 )(I BA - Ib,2) + Cii(0) + C 22 (0) - C 12 (0) - C 21 (0) 



Cab = 



(h/eV)(I A ,i + Ia,2)(Ib,i + Ib, 2 ) + Cn(0) + C 22 (0) + C 12 {0) + C 2 i(0) ' 

I 



(5.6) 



(5.7) 



(5.8) 



The correlator (|5.8|) differs from the expression (|5.5|) in the tunneling regime by the extra terms containing 
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the product of average currents. It still contains only 
low- frequency quantities, so it can be evaluated with- 
out requiring time-resolved detection. In Appendix IC.2I 
we show that the correlator l|5.8|l produces again the ex- 
pected result (13.11(1 for the concurrence, now without the 
requirement t\ , T2 <C 1 . 



C. Full counting statistics 

The relation between entanglement and current fluc- 
tuations has been extended to higher than second order 
correlators by Faoro, Taddei, and Fazio 68] and by Di 
Lorenzo and Nazarov |5J|. This is the problem of full 
counting statistics. A complication in the analysis is 
that electrons do not traverse a quantum conductor as 
separate and independent particles unless the transmis- 
sion probability is vanishingly small. There is no way to 
determine a priori if and how the electron flow can be 
partitioned into elementary events involving one, two, or 
larger number of electrons. In Ref. using the spin 
entangler of Fig. |Ht, it was shown that the elementary 
events of electron transfer which contribute to the full 
counting statistics can be viewed as single-electron and 
two-electron events. Two electrons are always transferred 
as the spin singlet (|2.1[1 . The fraction of electrons trans- 
ferred as singlets depends on the transmission probabil- 
ity. While all electrons come in singlets for a perfectly 
transmitting channel, no singlets are transferred in the 
limit of vanishing transmission probability. 



VI. LOSS OF ENTANGLEMENT BY DEPHASING 



In Sec. IIII.BI we discussed how thermal fluctuations 
degrade the entanglement of the electron-hole pair. The 
pure state at zero temperature becomes a mixed state 
as a result of the thermal fluctuations. Beyond a critical 
temperature T c of order eV/ks the entanglement produc- 
tion at the tunnel barrier vanishes identically (cf. Fig.[3J). 

If T <C eV/ks thermal fluctuations are irrelevant, but 
the electron-hole state may still have lost its purity (= 
have dephased or decohered) by other mechanisms. The 
two types of entanglement discussed in Sec. IIVI spin or 
orbital, dephase for different reasons. For spin entangle- 
ment, the spin-orbit interaction and the hyperfinc inter- 
action with nuclear spins are effective mechanisms of de- 
phasing. For orbital entanglement, the electromagnetic 
fluctuations caused by other charges lead to dephasing. 
In this section we concentrate on the latter case, specifi- 
cally in the edge channel geometry of Fig. ^jp. 

The dephasing by electromagnetic fluctuations of en- 
tangled electron-hole pairs in spatially separated edge 
channels has been studied in Refs. y, Wi |7(| . Here 
we follow Ref. |6|j. We model the effect of the electro- 
magnetic fluctuations phenomenologically by introducing 
phase factors in the scattering matrix and subsequently 
averaging over these phases. We refer to Ref. |7l| for an 



alternative phenomenological approach (using a dephas- 
ing voltage probe). A more microscopic treatment of the 
effect of dephasing on electron-hole entanglement (for ex- 
ample along the lines set out by Marquardt and Brudcr 
|72fl) has not yet been attempted. 

Introducing random phase shifts 4>i (ipi) accumulated 
in channel i at the left (right) of the tunnel barrier, the 
reflection and transmission matrices transform as 











ro, t 











to- 



(6.1) 



Pout 



By averaging over the phase shifts, with distribution 
P (61, fejipi, ^2), the pure state |* ut) [given by Eq. 
[4ll ] is converted into the mixed density matrix 

d(f>id4> 2 d^idtp2 P{(j)\, <fo,^i,^2)|*out)(*out|- 

(6.2) 

To simplify the average, the phase shifts at the left and 
the right of the tunnel barrier are assumed to be indepen- 
dent with identical distributions, so P(4>i, 4>%, ipi, V'a) — 
P(<t>i,<t>2)P{' l Pi,'4 > 2)- The complex dephasing parameter 
?7 is defined by 



r) = / d(j)d(j)' P '((f),(j)' ') e ^-^' '. 



(6.3) 



Analytical progress can be made in the case that the 
two transmission eigenvalues (eigenvalues of ttf) are iden- 
tical: Ti = T2 = t. In the absence of dephasing the elec- 
tron and hole then form a maximally entangled Bell pair, 
with concurrence C = 1. The matrix product r^Oyt^ , 



which appears in the expression 
rameterized by 



1 for |*out), is pa- 



roo-ytf = y/r(l-T) 



cos £ sin £ 
— sin £ cos £ 



(6.4) 



The angle £ governs the extent to which the tunnel bar- 
rier mixes the edge channels (no mixing for £ = 0, 7r/2, 
complete mixing for £ = 7r/4). 

Just as in the absence of dephasing (cf. App. fAJ the 
entanglement production is determined entirely by the 
projection of the density matrix p out on a single particle 
left and right of the tunnel barrier. The projected density 
matrix pn has elements 



/ 2cos 2 £ 



P11 



■q* sin 2£ 
— T] sin 2£ 
\2|r/| 2 cos 2 £ 



77 sin 2£ 
2 sin 2 £ 
-2t7 2 sin 2 £ 
rj sin 2£ 



-■q* sin2£ 
-2?7* 2 sin 2 £ 

2 sin 2 £ 
—rj* sin 2£ 



2 cos 2 A 



2\rj 
77* sin 2£ 
— rj sin 2£ 
2 cos 2 £ j 
(6.5) 

The weight W\\ = 2r(l — r) of the projection is indepen- 
dent of dephasing. 

From the mixed density matrix l|6.5|) we calculate the 
concurrence using Wootters formula (|2.19|l . with the re- 
sult 



C 



ix{ o,-i(i-H 2 ) 

Wl6M 2 + 2(l- 



2 (1 + cos4£)}. (6.6) 
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FIG. 11 The Bell inequality is violated for dephasing param- 
eters 1 7] | > rj' c , while entanglement is preserved for \r)\ > rj c . 
The dashed and dotted curves are computed from Eqs. I|ri.8[l 
and 16.101 1. respectively. The shaded region between the two 
curves indicates dephasing and mixing parameters for which 
there is entanglement without violation of the Bell inequality. 
Figure from Ref. |69l . 



Notice that C = \rj\ 2 for £ = 0. The entanglement pro- 
duction rate follows from 



£part / i 



dot 



(6.7) 



with the function T defined in Eq. 1)2. 12[) . 

Instead of a critical temperature we now have a criti- 
cal dephasing parameter rj Cl such that the entanglement 
production vanishes identically for |?7| < r\ c . (Smaller \q\ 
means stronger dephasing.) From Eq. l|6.6|l one finds the 
expression 



5-cos4£-2V6-2cos4£ 

^ = V r^oiii ' (6 - 8) 

plotted in Fig. (dashed curve). Unlike the critical 
temperature (j8.12|l . the critical dephasing parameter is 
independent of the tunnel probability r. It does depend 
sensitively on the mixing parameter £, increasing from 
to V2 — 1 as £ increases from to 7r/4. 

As discussed in Sec. III.Dl for a mixed state there is no 
one-to-one relation between the Bell parameter B max and 
the concurrence C. Application of the general formula 
(12.221) to the density matrix l|6.5[) gives 



Bmax = V / 2(l + h| 2 ) 2 + 2(l-|7 7 |2)2 C Os4e (6.9) 

The Bell inequality B max < 2 is violated for |?/| > r)' c , 
with 



, / - 1 + cos 4£ + V2 - 2 cos 4£ 
Vc = \ 



1 + cos 4£ 



(6.10) 



FIG. 12 Relation between the maximal violation S max of 
the Bell inequality and the concurrence C, calculated from 
Eqs. 133 and (EU for mixing parameters £ = (triangles, 
no mixing) and £ = 7r/4 (squares, complete mixing). The 
dephasing parameter |?7| decreases from 1 (upper right corner, 
no dephasing) to (lower left, complete dephasing) with steps 
of 0.05. The dotted line is the relation between £> max and C 
for a pure state, which is also the largest possible value of 
Bmax for given C. Figure from Ref. |o9l| . 



increasing from to V2 as £ increases from to 7r/4 
(dotted curve in Fig. Ill|) . 

In Fig. 1121 we compare B max and C for £ = (no mix- 
ing) and £ = 7r/4 (complete mixing). For £ = the 
same relation (I2.15|l between B max and C holds as for 
pure states (dotted curve). Violation of the Bell-CHSH 
inequality (|2.6(l is then equivalent to entanglement. For 
£ 7^ there exist entangled states (C > 0) without viola- 
tion of the Bell inequality (B max < 2). This is indicated 
by the shaded region in Fig. ^2 For £ ^ violation of 
the Bell inequality is thus a sufficient but not a necessary 
condition for entanglement. 



VII. QUANTUM ENTANGLEMENT PUMP 

The quantum electron pump is a device that trans- 
fers electrons phase coherently between two reservoirs at 
the same voltage, by means of a slowly oscillating volt- 
age on a gate electrode Special pump cycles exist 
that transfer the charge in a quantized fashion, one e per 
cycle. Samuelsson and Biittiker |74| asked the question 
whether a similar device could be used to create entan- 
gled electron-hole pairs in a controlled manner, clocked 
by the gate voltage, rather than stochastically as in the 
original proposal @] • (See Fig. ED for a schematic illus- 
tration.) Such a quantum entanglement pump could be a 
building block of quantum computing designs using bal- 
listic flying qubits in nanowires or in quantum Hall edge 
channels [zllzllz3 
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FIG. 13 Production of entangled electron-hole pairs in a 
narrow ballistic conductor by a quantum pump. The left 
and right ends of the conductor are at the same potential, 
while the potential on the gate electrodes at the center is pe- 
riodically modulated. Such a device produces spatially sep- 
arated electron-hole pairs (black and white circles), differing 
in energy by a multiple of the pump frequency u). For spin- 
independent scattering the electron (e) and hole (h) produced 
during a given cycle have the same spin f, 1, so that their wave 
function is that of a Bell pair, oc | TeTh) + | ieih)- The op- 
timal quantum entanglement pump produces, on average, 1 
Bell pair every 2 cycles. Figure from Ref. [7^1 . 



We already noticed in Sec. ITlLBl [cf. Eq. lpT5)l ] that 
the rate of entanglement production is closely related to 
the charge noise, to the extent that a noiseless scatterer 
produces no entanglement. By maximizing the charge 
noise with spin-independent scattering we calculated |78j 
that a pump can at most produce, on average, 1 Bell 
pair every 2 cycles. (A similar conclusion was reached 
in Ref. |z3.) A deterministic spin-entangler jgjj, being 
the analogue of a quantized charge pump, would have an 
entanglement production of 1 Bell pair per cycle, so the 
optimal entanglement pump has one half the efficiency of 
a deterministic entangler. 



VIII. TELEPORTATION BY ELECTRON-HOLE 
ANNIHILATION 

Teleportation is the disembodied transport of a quan- 
tum mechanical state between two locations that are only 
coupled by classical (incoherent) communication |81| . 
What is required is that the two locations share a previ- 
ously entangled state. Teleportation has the remarkable 
feature that the teleported state need not be known. It 
could even be undefined as a single-particle state, which 
happens if the teleported particle is entangled with an- 
other particle that stays behind. Teleportation then leads 
to "entanglement swapping" [82L l83| : Pre-existing en- 
tanglement is exchanged for entanglement between two 
parties that have never met. 

Experiments with photons [84| have demonstrated that 
teleportation can be realized in p ractice. Only linear op- 
tical elements are needed [85| 86] , if one is satisfied with a 
success probability < 1. Such non-deterministic telepor- 
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FIG. 14 Schematic description of teleportation by electron- 
hole annihilation. A voltage V applied over a tunnel barrier 
produces pairs of entangled electron-hole pairs in the Fermi 
sea. One such pair (e_L,/i_L) is shown at the left. For a sim- 
plified description we assume spin entanglement in the state 
(I TT) + I ll))/V2, where the first arrow refers to the electron 
spin and the second arrow to the hole spin. A second electron 
en is in an unknown state ct\ f) + /3| [). The electron en can 
annihilate with the hole by tunneling through the barrier 
at the center. If it happens, the state of eL collapses to the 
state of en- (Notice that | f ) annihilates with | f) and j J.) an- 
nihilates with ||), so eL inherits the coefficients a and of en 
after its annihilation.) The diagram shows a second entangler 
at the right, to perform two-way teleportation (from e_R to e_L 
and from fti to ha). This leads to entanglement swapping: 
e_L and ha become entangled after the annihilation of Kl and 
e_R. Figure from Ref. |88|. 



tation plays an essential role in proposals for a quantum 
computer based entirely on linear optics |25j . A cen- 
tral requirement for nontrivial logical operations is that 
the linear elements (beam splitters, phase shifters) are 
supplemented by single-photon sources and single-photon 
detectors, which effectively introduce nonlinearities. 

Teleportation of electrons has not yet been realized. 
The analogue of teleportation by linear optics would be 
teleportation of free electrons, that is to say, teleporta- 
tion using only single-particle Hamiltonians. A direct 
translation of existing linear optics protocols would re- 
quire single-electron sources and single-electron detectors 
|87|. Such devices exist, but not for free electrons — they 
are all based on the Coulomb interaction in quantum 
dots. 

In this section we review, following Ref. 88], an alter- 
native single-particle teleportation mechanism using en- 
tangled electron-hole pairs. The key observation is that 
the annihilation of an electron-hole pair in the Fermi sea 
teleports these quasiparticles to a distant location, if en- 
tanglement was established beforehand. This two-way 
teleportation scheme is explained in Fig. 

The analysis is simplest for the entangled state 

(|T>e|t>A + U>eU>A)/vl 

The subscripts e and h refer, respectively, to electron and 
hole at two distant locations. The particle to be tele- 
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FIG. 15 Proposed realization of the teleportation scheme of Fig. 1141 using edge channels in the quantum Hall effect. The 
thick black lines indicate the boundaries of a two-dimensional electron gas, connected by Ohmic contacts (black rectangles) to 
a voltage source V or to ground. A strong perpendicular magnetic field ensures that the transport in an energy range eV above 
the Fermi level takes place in two edge channels, extended along a pair of equipotentials (thin solid and dashed lines, with 
arrows that give the direction of propagation) . These edge channels realize the two-channel conductors of Fig. 1141 with the 
Landau level index n — 1, 2 playing the role of the spin index f , J.. Solid lines signify predominantly filled edge channels with 
hole excitations (open circles), while dashed lines signify predominantly empty edge channels with electron excitations (black 
dots). The beam splitters of Fig. 1141 are formed by split gate electrodes (shaded rectangles), through which the edge channels 
may tunnel (dashed arrows, scattering matrices Sl, Sr, So). The annihilation of the electron-hole excitation at the central beam 
splitter is detected through the currents I a and Ib- Entanglement swapping resulting from two-way teleportation is detected 
by the violation of a Bell inequality. This requires two gate electrodes to locally mix the edge channels (scattering matrices Ul , 
Ur) and two pair of contacts 1, 2 to separately measure the current in each transmitted and reflected edge channel. Notice that 
there are no phase-coherent paths connecting the left and right ends of the conductor (because of the intervening dephasing 
contacts A and B), so a demonstration of entanglement between the two ends is indeed a demonstration of teleportation. Figure 
from Ref. 0. 



ported is another electron, in the state a\ 1) e > + j3\ |) e ' 
(with |a| 2 + |/3| 2 = 1). The second electron e' may tunnel 
into the empty state representing the hole h, but only if 
the spins match. If r denotes the tunneling amplitude, 
then this happens with probability !M 2 |t| 2 -|-^|/3| 2 |'?"| 2 — 
\\t\ 2 <C 1. The resulting annihilation of the two quasi- 
particle excitations collapses the combined state 

HT)e'+/3|l}eO(IT}e|T)h + |l}e|4)h)A/2 

to the state a\ T) e + P\ l)e, so the state of the second 
electron e' is teleported to the first electron e at a distant 
location. 

The usual limitations [8 If of teleportation apply. Since 
tunneling is an unpredictable stochastic event, it has to 
be detected and communicated (by classical means) to 
the distant location. There is therefore no instantaneous 
transfer of information. And since the electron has to be 
annihilated in order to be teleported, its state can not 
be copied. Teleportation by electron-hole annihilation 
thus presents a rather dramatic demonstration of the no- 
cloning theorem of quantum mechanics [8^, |9(j ■ 

A major obstacle to teleportation in the solid state is 
the requirement of fast time-resolved detection. To cir- 
cumvent this difficulty Ref. |8S| identified a low-frequency 
noise correlator that demonstrates the entanglement 



swapping resulting from two-way teleportation. Two-way 
teleportation means that upon annihilation the electron 
and hole are teleported to opposite ends of the system. 
The noise correlator measures the degree of entanglement 
at the two ends. This demonstrates teleportation if the 
two ends are not connected by any phase-coherent path. 

In Fig. 1151 a particular implementation is illustrated 
using edge channels in the quantum Hall effect regime. 
To detect the entanglement swapping one needs to mea- 
sure a fourth order cumulant of fluctuations of tunneling 
current. Typically, only the second order cumulant is 
measured in noise experiments, but higher order cumu- 
lants are now becoming accessible as well 0, |9^ • 



IX. THREE-QUBIT ENTANGLEMENT 

So far we have always considered bipartite entangle- 
ment. The interaction-free entanglement of three qubits 
(= tripartite entanglement) was investigated in Ref. |93| . 

The proposed three-qubit ent angler is sketched 
schematically in Fig.^J As in the original three-photon 
entangler of Zeilinger et al. |94| . three-qubit entangle- 
ment is created out of two entangled electron-hole pairs. 
Fig.ll7lshows a physical realization of Fig. ^1 using edge 
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FIG. 16 Schematic description of the creation of three- 
qubit entanglement out of two entangled electron-hole pairs 
in the Fermi sea. The left and right entangler consist of 
a tunnel barrier over which a voltage V is applied. For a 
simplified description we assume spin entanglement in the 
state (| Tfcte) + | |fcie))/v / 2 I where the subscripts e,h re- 
fer to electron and hole spin. The two electrons meet at a 
polarizing beam splitter, which fully transmits the up-spin 
and fully reflects the down-spin. If the outgoing ports A, 
B contain one electron each, then they must both have the 
same spin. The corresponding outgoing state has the form 
(I T/itftTeTe) + | |ft|h|e|e))/V2- Since the two electrons at 
A,B are constrained to have the same spin, this four-particle 
GHZ state represents three independent logical qubits. Figure 
from Ref. 1931. 



channels in the quantum Hall effect. 

The irreducible tripartite entanglement is quantified 
by the tangle r of Coffman, Kundu, and Wootters 
[95j, which is the three-qubit analogue of the concur- 
rence. The tangle is unity for the maximally entangled 
Greenberger-Horne-Zeilinger (GHZ) state | ttt) + I ill) 
and vanishes if one qubit is disentangled from the other 
two. Unlike the concurrence (|3.11|) . the tangle depends 
not only on the transmission eigenvalues of the point con- 
tact entanglers, but also on the eigenvectors. 

As in the bipartite case, current correlators can be used 
to test for the maximal violation of a tripartite general- 
ization of the Bell inequality. Two tripartite inequali- 
ties are known, one due to Mermin and the other 
to Svetlichny [93. While there exists a one-to-one rela- 
tion (|2.15|l between concurrence and Bell inequality for 
any pure state of two qubits, no such relation exists for 
r. A numerical investigation |98j| has found a simple set 
of upper and lower bounds for r. Since these bounds 
become tighter and tighter as the state approaches the 
GHZ state, they should be of practical use. 

The bipartite electron-hole entangler [3j discussed in 
previous sections is capable of producing the most gen- 
eral two-qubit entangled pure state, by suitably choosing 
the scatte ring matrix of the tunnel barrier. The tripartite 
entangler [93 discussed here, in contrast, is limited to the 
production of a three-parameter subset of the most gen- 
eral five-parameter family of three-qubit entangled pure 
states [93. This subset is characterized by the prop- 
erty that tracing over the third qubit results in a mixed 



two-qubit state which is not entangled. The origin of this 
restriction is that the three-qubit state is constructed out 
of two separate entangled electron-hole pairs. 

The two-qubit entangler can produce maximally en- 
tangled Bell pairs as well as partially entangled states, as 
quantified by the concurrence. Similarly, the three-qubit 
entangler can produce maximally entangled GHZ states 
as well as states that have a smaller degree of tripartite 
entanglement, as quantified by the tangle. However, in 
the three-qubit case there is a second class of states that 
are irreducibly entangled and can not be ob tained from 
the GHZ state by any local operation [lOOf . These so- 
called W-states are not accessible by the method of Fig. 
IT51 It would be interesting to see if there exists an al- 
ternative interaction-free method to extract the W-state 
out of the Fermi sea, or whether this is impossible as a 
matter of principle. 



X. THE EXPERIMENTAL CHALLENGE 

Tunneling experiments in metals and semiconductors 
have been carried out in the laboratory since the late 
1950's, following the pioneering work of Esaki and Gi- 
aever \m\ . The prediction of Ref. is that at low tem- 
peratures these experiments produce electron-hole pairs 
in the spin-entangled state l|2.1|) . So we may well have 
been producing Bell pairs in the solid state even earlier 
than in optics! The rea son t hat the Bell inequality has 
been violated in optics |l02j| but not yet in electronics 
seems, therefore, to be a difficulty of detection rather 
than production. It is indeed much easier to measure 
photon polarization than electron spin. 

Some of the implementations discussed in Sec. IIVI do 
not require spin-resolved detection, but since the entan- 
glement then involves orbital degrees of freedom it is 
much more sensitive to decoherence. Attempts to ob- 
serve two-particle interference effects of charge currents 
(for example, in the Mach-Zehnder geometry of Ref. |6t| ) 
have so far been unsuccesful, presumably because the 
electr on ch arge is too strongly coupled to the environ- 
ment ,1031 . Average spin currents can now be reliably 
detected |l04j . It seems likely that correlators of spin cur- 
rents will become accessible experimentally in the near 
future, permitting a test for the predicted electron-hole 
entanglement in the Fermi sea. 
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FIG. 17 Proposed realization of the three-qubit entangler, using edge channels in the quantum Hall effect. The left and 
right point contacts (scattering matrices Sl, Sr) each produce entangled electron-hole pairs in the Fermi sea. They partially 
transmit and reflect both edge channels, analogously to beam splitters in optics. The central point contact is the analogue of a 
polarizing beam splitter: It fully transmits the inner edge channel and fully reflects the outer one. Three-qubit entanglement 
results if there is one excitation at each of the four edges L, R, A, B. The two electron excitations at A and B then have 
the same channel index, so they constitute a single qubit. This qubit forms a three-qubit entangled state with the two hole 
excitations at L and R. Figure from Ref. |93l| . 



APPENDIX A: Entanglement production for 
spin-dependent scattering 

In this appendix we review the calculation leading to 
Eqs. (|3~TT1 and jXTJl , following Ref. Q. 

At zero temperature the incident state has the form 

l*in)= LI 4,i(£)4, 2 vE0|O>- (Al) 

E F <E<E F +eV 

The fermion creation operator a\ n (E) (with n = 1, 2 =1, 
J.) excites the n-th channel incident from the left at en- 
ergy E. Similarly, a, Rn (E) excites a channel incident 
from the right. The product over energies refers to a dis- 
cretization with interval 8E = h/tdet- The total number 
of energy intervals that contribute is eV/SE — eVt^ct/h. 

It is convenient to collect the creation operators in two 
vectors a L , a R and to use a matrix notation, 

i^-n(:|) (*?!!)(:{>■ w 

The input-output relation (|3 . H|> can then be written as 

fc) - (; .''*) (s) ■ 

The 4x4 unitary scattering matrix S has 2x2 submatri- 
ces r,r',t,t' that describe reflection and transmission of 
states incident from the left or from the right. The trans- 
mission eigenvalues t\,t% S [0,1] are the eigenvalues of 
tH = 1 — r^r (which are the same as the eigenvalues of 
t'H' = 1 — r''r' because of unitarity) . Substitution of Eq. 



(|A3|I into Eq. i|A2(l gives the outgoing state 



|*out) = Y[{b{ra y t T b R + [ra y r T } 12 b[ A bl 2 

E 

+ [ta y t T } 12 b RA b R2 )\0). (A4) 



The superscript T indicates the transpose of a matrix. 

As discussed in Sec. III. El to calculate the accessible 
entanglement we need to account for the constraint of 
local particle number conservation in each energy interval 
SE. 7 The three terms at the right-hand-side of Eq. (|A4|) 
represent different local particle numbers p and q at the 
left and right of the tunnel barrier, so the entanglement 
entropy of each term may be added to obtain the total 
entanglement entropy £ par t of the outgoing state [cf. Eq. 
(I2.25j . The second term (p = 2, q = 0) and the third 
term (p = 0, q = 2) are separable states. Only the first 



7 One could relax this constraint, requiring only that the total 
particle number, summed over all energies, is conserved at the left 
and at the right of the tunnel barrier. The constraint used here, 
of particle conservation in each energy interval <5E, is appropriate 
for an observer who uses only elastic scattering to detect the 
entanglement production. 
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term (p = q = 1) gives a nonzero contribution, 8 
eVt dct 



part 



(I*)), 



(A5a) 

|$) = 6+7^10), 1 = w- l / 2 ra y t T , (A5b) 
w = Tr (r<7 y t T )(ra y t T )1 = n + r 2 - 2nT 2 .(A5c) 

The factor eW<jet//i counts the number of independently 
contributing energy intervals. 

The degree of entanglement of the two-particle pure 
state |$) is quantified by the concurrence. Substituting 
Eq. (|A5|> into Eq. (|2.13|) one arrives at the result 
The result (I3.10|) for the entanglement entropy follows 
from Eq. (t2~TT)l . 



APPENDIX B: Entanglement production at finite 
temperature 

In this appendix we calculate the free-fermion entan- 
glement produced at a tunnel barrier for the case that 
the electron sources are at finite temperature, general- 
izing the zero-temperature calculation of Ref. This 
calculation not been previously published. 

We refer to the beam splitter geometry of Fig. f|] The 
incoming state p ln produced by the sources at an elevated 
temperature and/or voltage is scattered via a beam split- 
ter into an outgoing state p ut, which is absorbed by de- 
tectors at zero temperature and voltage. (The detectors 
should be kept at a lower temperature than the sources, 
in order to be able to distinguish incoming from outgoing 
states.) For electron sources in local thermal equilibrium, 
the density matrix has the form 



nn 



(l-f n (E))\0)(0\+f n (E)aUE)\0)(0\a n (E) 



( B1 ) 

Here a n (E) is the fermion annihilation operator for an 
incoming channel n at energy E and f n {E) is the Fermi- 
Dirac distribution. (Different channels may have a differ- 
ent temperature and/or voltage, which is why we speak 
of local rather than global thermal equilibrium.) The 
state |0) is the vacuum state (such that a„|0) =0). 

8 In the case of spin-independent scattering, when 7 = 
2~ 1 / 2 a y , the state |$) simplifies to |f>) = -i2~ 1 / 2 (b , L 1 6j J 2 - 

bl i2 b ] R1 )\V) = -i2- 1 / 2 (| T)e| |>e " I i)e| T)e), which is the 

spin singlet state of Eq. 12.11 . The electron- hole Bell pair 



The annihilation operators b n (E) for the outgoing 
channels are related to those for the incoming channels 
by the unitary scattering matrix of the tunnel barrier, 



b n (E) 



(B2) 



/ J Snm{E)a m (E). 

Substitution of Eq. (|B2|1 into Eq. IjBlfl gives the density 
matrix of the outgoing state in a given energy interval 
SE, 



Pout 



n 



(i-/„)|o)<o|+/„4 



"n — y ] bln^n 



(B3) 

We consider the case of a single spin-degenerate mode, 
so n = 1,2 labels spin up and down from source A and 
n = 3,4 labels spin up and down from source B. We 
denote fx = f% = Ja and /a = /4 = /s- The sources are 
taken at the same temperature T with a voltage differ- 
ence V, 



Ia{E) 



1 



exp[(£ -E F - eV)/k B T] + 1 ' 

ME) = eM{E-E l F )lk B T] + l- (B4) 

The Fermi energy at zero voltage is denoted Ep. 

To determine the entanglement production, con- 
strained by particle conservation, we need to project the 
density matrix p out onto sectors N pq of Fock space with a 
well-defined particle number p,q £ {0, 1, 2} at detector A 
and B in each energy interval SE. In accordance with Eq. 
(|2.25|l , the entanglement production £ par t is the weighted 
sum of the entanglement in the individual sectors. The 
only sector with nonzero entanglement is A/"n, with one 
particle at detector A and one particle at detector B. In 
all other sectors (with pq — 0, mod 2) the projection p pq 
is separable in the two degrees of freedom at detectors A 
and B. 

The projected density matrix pxi is constructed from 
Eq. HB3fl by summing over the six nonequivalent ways 
to partition the indices 1,2,3,4 into two filled channels 
nx, ri2 and two empty channels 713,714, 

2 _1/2 (| T)fc| T)e +| l)h\l)e) of Eq. 03 is obtained by a canon- 
ical transformation from electron to hole degrees of freedom at 
the left of the tunnel barrier. 



wxxpn = fnJnA^ - fn,)(l ~ fn^ A M(n 1 ,n 2 )b B \0)(0\b B MHn 1 ,n 2 )b A , (B5) 

partitions 

wxx = fnJn 2 (l-fn s )(l-fn A )TrM(nx,n 2 )M^nx,n 2 ), (B6) 

partitions 

£>3m <S , 3n 2 j (B7) 



M{ ni ,n 2 ) = S q ln Aa 
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We have defined the vectors bA = (£>i, 62), &b = (63, 64) of 
annihilation operators at detectors A and B, respectively. 

The general solution simplifies considerably if the scat- 
tering matrix S is spin- independent. It then has the block 
diagonal form 



the entire energy dependence of <5£ par t lies in the weight 
factor wn. Integration over energy gives the total entan- 
glement production 



eVt dc 



part 



■T 



3 i(</>+0')l 



2V T 



(B8) 



cosh -i 



l-±(3£-l) 2 
- u^ 1 sinhu 



(1 - shfiT 



(B18) 



with I2 the 2x2 unit matrix. We have parameterized 
the transmission and reflection amplitudes in terms of 4 
phases <j>,<f/,ip,ip', and the transmission probability r S 
[0,1]. 

The projected density matrix l|B5|) becomes 

w lxPll = Ai (| TTXTT I + 1 11X11 1) + A 2 |*Boii)(*Bcii| 



^3(l*7)(*n 



(B9) 

i«u = 4/ a (1-/ a )/ b (1-/ s ) + 2(/ a -/ b ) 2 t(1-t), 

(BIO) 

Ai = /a(1-/a)/b(1-/b), (Bll) 
A 2 = 2[/ 2 (l-/ B ) 2 +/ 2 (l-/ A ) 2 ]r(l-r), (B12) 
A3 = /. 4 (1 - /a)/b(1 - - 2r(l - r)]. (B13) 

We denote the spin channels by t,l, so that, for exam- 
ple, b\b\\0) = I H). In Eq. lp))l. the first term oc Ai 
contributes two separable states, the second term oc A2 
contributes the maximally entangled Bell state Q2.1JI . and 
the third term oc A3 contributes a mixture of the partially 
entangled state 

|* 7 ) = [1 - 2r(l - r)]-^[(l - r)| U) + r| |t>] (B14) 

and its spin-inverse. 

With some more algebra the density matrix (|B9|) re- 
duces to the Werner state l|2.20|l . 



This function is plotted in Fig. 

The maximal entanglement production rate of eV/2h 
bits per second is reached for T = 0, r = 1/2. At 
nonzero temperature the entanglement production de- 
creases, vanishing oc (T c — T) 2 hi(T c — T) at the critical 
temperature T c determined by Eq. Ij3.12|l . 



APPENDIX C: Bell inequality with noise correlators 
1. Tunneling regime 

Following Ref. we calculate the maximal violation 
of the CHSH inequality l|2.6|l that follows using the ex- 
pression l|5.5|l for the spin-spin correlator C a & in terms of 
spin-resolved noise correlators. 

At low temperatures (fc^T <C eV) the zero-freque ncy 
noise correlator dj (0) has the general expression j^J, LlJ 

M 



C i3 = -(e 3 V/h)\(rtf) 



ij 



(CI) 



PU = i(l-£)l4 + £l*Bcn)(*Boll| 

£, = (2/ Wll )(fA-fB) 2 r(l-r) 
= 1 



1 



1 + 2t(1 - t) sinlr it 



(B15) 



(B16) 



We have abbreviated u = eV/2ksT. Notice that the 
energy dependence introduced by the Fermi-Dirac dis- 
tribution drops out of the final expression for £. The 
concurrence of the Werner state is given by Eq. H2.21JI . 
Substitution into Eqs. Ij2.18|l and (|2.25(l gives the entan- 
glement production 5£ vax t in the energy interval 5E un- 
der consideration, 



in terms of the 2x2 reflection and transmission matrices 
r, t of the tunnel barrier [cf. Eq. (|A3|) ]. It is convenient to 
parameterize these matrices in terms of the transmission 
eigenvalues ri , t% , by 

(C2) 

with 2x2 unitary matrices Ul, Ur, Uq- [In the presence 
of time-reversal symmetry one has Uq = lf£ ' , but this 
constraint is irrelevant because anyway Uq drops out of 
Eq. KH1 .1 

The spin-spin correlator C a b depends on the directions 
a and b along which the spin is measured at the left and 
right of the tunnel barrier. Combination of Eqs. I|5.5I) 
and ijClfl gives 



Cab = 



(B17) 

with 9(x) the unit step function \9(x) = if x < 0, 
6{x) = 1 if x > 0]. 

We assume that the transmission probability is approx- 
imately energy independent on the scale of fcsT 1 . Then 



Tr (a • cr)rtf (b ■ cr)tr^ 
Tr rtHr^ 



(C3) 



Equivalently, one can measure the spins along the fixed 
z-axis after having rotated them by unitary matrices U a 
and Ub, such that 



U a a z U a = a-<r, Ula z U b = b-a. (C4) 
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In the parameterization (|C2|) one finds 

C ab = (l-2|?7 a41 | 2 )(l-2|c7 b!ll | 2 ) 

+ 4nCRcU aA1 UZ M U: A2 U bA2 , (C5) 
U a = U a U L , U b = U b U R , (C6) 

K=l+ „ ( V T2) l r. (C7) 

n(i -n) + t 2 (i -r 2 ) 

The quantity C in Eq. I|C5|I is the concurrence (|ihll|) of 
the electron-hole pair. 

By varying the unit vectors a, b, a', b' (or, equivalently, 
the unitary matrices U a , U b , U a ', U b >) in the CHSH in- 
equality (|2.ti|) one finds the maximum violation 

B max = 2^1 + k 2 C 2 m 2y/l +C 2 if n, r 2 < 1. (C8) 

This confirms the expected relation ({2.15(1 between the 
concurrence and the Bell parameter in the tunneling 
regime t\,t% <C 1. 

2. Beyond the tunneling regime 

Following Ref. [H^, we calculate the maximal viola- 
tion of the CHSH inequality (|2.6|) that follows using the 
expression (|5.8|l for the spin-spin correlator C ab . This ex- 
pression contains both the noise correlators and the mean 
currents and does not require the tunneling approxima- 
tion. Referring to the beam splitter geometry of Fig. 0] 
only the current from the beam splitter to the detectors 
should be included in the mean current, not the current 
from the sources to the beam splitter 9 . 

The mean currents into the detectors are given by 

(lA,i) = ^(TTt)„, (j Bi ) = f^( tt t) ( C 9) 
h h 

Combination of Eqs. (|5.8p. (|Clfl . and (|C9|) gives 

_ [Tr (a ■ er)rrt] [Tr (b ■ tr)tt t ] - Tr (o • <r)rt T (6 ■ cr)tr^ 
ab ~ (Tr rr t)(Trttt) _ Tr rtHr^ 

(CIO) 

Substitution of Eqs. I|C2|I and (|C4|) leads to 

C ab = -(l-2|L7 aill | 2 )(l-2|c\ xl | 2 ) 

- 4CReU a>n U^ n U* >12 U bA2 . (Cll) 

Apart from an overall minus sign, Eq. (|C11|) is the same 
as Eq. I(C5() without the spurious factor n multiplying the 



This distinction between outgoing and incoming particles is ir- 
relevant in the tunneling regime at zero temperature, because 
the non-fluctuating incoming current does not contribute to the 
noise correlator dj(0). In an electronic circuit the outgoing cur- 
rent can be measured separately from the incoming current in a 
four-terminal geometry with grounded detectors at zero temper- 
ature. 



concurrence. The maximal violation of the Bell inequal- 
ity becomes 

£ max = 2^1 + C 2 . (C12) 

The relation (|2.15(1 between concurrence and Bell param- 
eter is now satisfied beyond the tunneling regime. 
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